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Abstract 

Strong electrostatic turbulence in magnetically confined plasmas is charac- 
terized by trapping or eddying of particle trajectories produced by the E x B 
stochastic drift. Trapping is shown to produce strong effects on test particles 
and on test modes. It determines non-standard statistics of trajectories: non- 
Gaussian distribution, memory effects and coherence. Trapped trajectories 
form quasi-coherent structure. Trajectory trapping has strong nonlinear ef- 
fects on the test modes on turbulent plasmas. We determine the growth rate 
of drift modes as function of the statistical characteristics of the background 
turbulence. We show that trapping provides the physical mechanism for the 
inverse cascade observed in drift turbulence and for the zonal flow generation. 
Keywords: plasma turbulence, statistical approaches, test particle transport, 
Lagrangian methods 

1 Introduction 

A component of particle motion in magnetized plasmas is the stochastic electric 
drift produced by the electric field of the turbulence and by the confining magnetic 
field. This drift determines a trapping effect or eddy motion in turbulence with slow 
time variation [1]. Typical particle trajectories show sequences of trapping events 
(trajectory winding on almost closed paths) and long jumps. Numerical simulations 
have shown that the trapping process completely changes the statistical properties 
of the trajectories. Particle motion in a stochastic potential was extensively studied 
[2]-[4], but the process of trapping was not described until recently. 

New statistical methods were developed [5], [6] that permitted to determine the 
effects of trapping. These are semi-analytical methods based on a set of deterministic 
trajectories obtained from the Eulerian correlation of the stochastic velocity. It was 
shown that trapping determines memory effects, quasi- coherent behavior and non- 
Gaussian distribution [6]. The trapped trajectories have quasi-coherent behavior and 
they form structures similar to fluid vortices. The diffusion coefficients decrease due 
to trapping and their scaling in the parameters of the stochastic field is modified. We 
have shown that anomalous diffusion appears due to collisions and average flows. A 
review of the effects of trapping on test particle statistics and on turbulent transport 
is presented in the first part of this paper. 
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The effects of trajectory trapping on tlie nonlinear dynamics of tlie test modes 
for tlie drift turbulence are presented in the second part of the paper. The semi- 
analytical methods developed for test particles are extended to test mode evolution 
in a turbulent magnetized plasmas. Test modes are usually studied for modelling 
wave- wave interaction in turbulent plasmas [7]. A different perspective is devel- 
oped here by considering test modes on turbulent plasmas. They are described by 
nonlinear equations with the advection term containing the stochastic E x B drift 
described by a stochastic field with known statistical characteristics. The growth 
rate of the test modes is determined as function of these statistical parameters. We 
develop a Lagrangian approach of the type of that introduced by Dupree [8], [9]. 
The difference is that in Dupree's method the stochastic trapping of trajectory was 
neglected and consequently the results can be applied to quasilinear turbulence. Our 
method takes into account the trapping and the non-standard statistics of trajecto- 
ries that it yields and thus it is able to describe the nonlinear effects appearing in 
strong turbulence. 

The paper is organized as follows. The test particle model is presented in Sec- 
tion 2. Section 3 contains and a short description of the statistical methods. The 
nonlinear effects of trajectory trapping on test particle statistics and transport are 
presented in Section 4. The general physical explanation for the anomalous diffusion 
regimes appearing in the presence of trajectory trapping and the formation of trajec- 
tory structures are discussed in this section. The problem of test modes in turbulent 
plasmas for the case of drift turbulence is presented in Section 5 where the growth 
rate and the frequency are determined as function of the statistical characteristics 
of the turbulence. The complex effects of trajectory trapping on the drift modes are 
analyzed in Section 6. The conclusions are summarized in Section 7. 



The test particle studies rely on known statistical characteristics of the stochastic 
field. They are determined from experimental studies or numerical simulations. The 
main aim of these studies is to determine the diffusion coefficients. The statistics 
of test particle trajectories provides the transport coefficients in turbulent plasmas 
without approaching the very complicated problem of self-consistent turbulence that 
explains the detailed mechanism of generation and saturation of the turbulent poten- 
tial. The possible diffusion regimes can be obtained by considering various models 
for the statistics of the stochastic field. 

We consider in slab geometry an electrostatic turbulence represented by an elec- 
trostatic potential 0^(x, t), where x = {xi,X2) are the Cartesian coordinates in the 
plane perpendicular to the confining magnetic field directed along z axis, B = -Be^. 
The test particle motion in the guiding center approximation is determined by 



where x(t) represent the trajectory of the particle guiding center, V is the gradient 
in the {xi,X2) plane and 0(x, t) = 0'^(x, The electrostatic potential 0(x, t) is 

considered to be a stationary and homogeneous Gaussian stochastic field, with zero 



2 Test particle model 
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average. It is completely determined by the two-point Eulerian correlation function 
(EC), E{x,t), defined by 

£;(x, t) = (0(x', t') 0(x' + X, + t)) . (2) 

The average (...) is the statistical average over the realizations of 0(x, t), or the 
space and time average over X' and P. This function evidences three parameters 
that characterize the (isotropic) stochastic field: the amplitude $ = •\/£'(0, 0), 
the correlation time Tc, which is the decay time of the Eulerian correlation and 
the correlation length Ac, which is the characteristic decay distance. These three 
parameters combine in a dimensionless Kubo number 

K = rjTfi (3) 

where tji = Xc/V is the time of flight of the particles over the correlation length 
and V — ^/Xcis the amphtude of the stochastic velocity. 
The diffusion coefficient is determined as (see [10]) 

A(i)= fdrLuir) (4) 
Jo 

where 

Lij(t; h) = {vi(O,0)vj{:>c(t),t)) (5) 

is the correlation of the Lagrangian velocity (LVC). It is obtained using the decor- 
relation trajectory method, a semi-analytical approach presented below. 

Equation (1) represents the nonlinear kernel of the test particle problem. The 
statistical methods will be presented for Eq. (1) for simplicity. They were developed 
to include complex models with other components of the motion (particle collisions, 
average flows, motion along the confining magnetic field, etc.). The effects of these 
components on the transport will be discussed in Section 4. 



3 The nested subensemble approach 

Test particle transport in magnetized plasmas in the nonlinear regime characterized 
by trajectory trapping was analytically studied only the last decade when the decor- 
relation trajectory method (DTM) [5] and the nested subensemble approach (NSA) 
[6] were developed. Trajectory trapping is essentially related to the invariance of 
the Lagrangian potential. Thus, a statistical method is adequate for the study of 
this process if it is compatible with the invariance of the potential. The NSA is 
the development of the DTM as a systematic expansion that obtains much more 
statistical information. 

The main idea in NSA is to study the stochastic equation (1) in subensembles of 
realizations of the stochastic field. First the whole set of realizations R is separated 
in subensembles (5*1), which contain all realizations with given values of the potential 
and of the velocity in the starting point of the trajectories x = 0, t = 0: 

(SI): 0(O,O) = (/)°, v(O,0) = v°. (6) 
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Then, each subensemble (SI) is separated in subensembles {S2) corresponding to 
fixed values of the second derivatives of the potential in x = 0, t = 



{S2): 0,^.(0,0) 



dxidxj 



x=O,t=0 
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where ij = 11,12,22. Continuing this procedure up to an order n, a system of 
nested subensembles is constructed. The stochastic (Eulerian) potential and veloc- 
ity in a subensemble are Gaussian fields but non-stationary and non-homogeneous, 
with space and time dependent averages and correlations. The correlations are zero 
in X = 0, t = and increase with the distance and time. The average potential and 
velocity performed in a subensemble depend on the parameters of that subensemble 
and of the subensembles that include it. They are determined by the Eulerian cor- 
relation of the potential (see [6] for details). The stochastic equation (??) is studied 
in each highest order subensemble (Sn). The average Eulerian velocity determines 
an average motion in each (Sn). Neglecting the fluctuations of the trajectories, the 
average trajectory in {Sn), X(t; Sn), is obtained from 



This approximation consists in neglecting the fluctuations of the trajectories in the 
subensemble {Sn). It is rather good because it is performed in the subensemble {Sn) 
where the trajectories are similar due to the fact that they are super-determined. 
Besides the necessary and sufficient initial condition x(0) = 0, they have supplemen- 
tary initial conditions determined by the definition (6-7) of the subensembles. The 
strongest condition is the inial potential 0(0, 0) = 0° that is a conserved quantity in 
the static case and determines comparable sizes of the trajectories in a subensem- 
ble. Moreover, the amplitude of the velocity fluctuations in {Sn), the source of the 
trajectory fluctuations, is zero in the starting point of the trajectories and reaches 
the value corresponding to the whole set of realizations only asymptotically. This 
reduces the differences between the trajectories in {Sn) and thus their fluctuations. 

The statistics of trajectories for the whole set of realizations (in particular the 
LVC) is obtained as weighted averages of these trajectories X(t; Sn). The weighting 
factor is the probability that a realization belongs to the subensemble {Sn); it is 
analytically determined. 

Essentially, this method reduces the problem of determining the statistical be- 
havior of the stochastic trajectories to the calculation of weighted averages of some 
smooth, deterministic trajectories determined from the EC of the stochastic poten- 
tial. This semi-analytical statistical approach (the nested subensemble method) is 
a systematic expansion that satisfles at each order n > 1 all statistical conditions 
required by the invariance of the Lagrangian potential in the static case. The order 
n — 1 corresponds to the decorrelation trajectory method introduced in [5] . In this 
case only the average potential is conserved. 

The nested subensemble method is quickly convergent. This is a consequence of 
the fact that the mixing of periodic trajectories, which characterizes this nonlinear 
stochastic process, is directly described at each order of our approach. The results 
obtained in flrst order (the decorrelation trajectory method ) for D{t) are practically 



dX{t; Sn) _ d^{X; Sn) 
It ~ dXj 



(8) 



A 



not modified in the second order [6]. Thus, the decorrelation trajectory method 
is a good approximation for determining diffusion coefficients. The second order 
nested subensemble method is important because it provides detailed statistical 
information on trajectories: the probability of the displacements and of the distance 
between neighboring trajectories in the whole ensemble of realizations and also in 
the subensembles (5*1). A high degree of coherence is so evidenced in the stochastic 
motion of trapped trajectories. 

4 Trapping effects on test particles 

4.1 Trajectory structures 

Detailed statistical information about particle trajectories was obtained using the 
nested subensemble method [6]. This method determines the statistics of the tra- 
jectories that start in points with given values of the potential. This permits to 
evidence the high degree of coherence of the trapped trajectories. 

The trapped trajectories correspond to large absolute values of the initial po- 
tential while the trajectories starting from points with the potential close to zero 
perform long displacements before decorrelation. These two types of trajectories 
have completely different statistical characteristics [6]. The trapped trajectories 
have a quasi-coherent behavior. Their average displacement, dispersion and proba- 
bility distribution function saturate in a time r^. The time evolution of the square 
distance between two trajectories is very slow showing that neighboring particles 
have a coherent motion for a long time, much longer than r^. They are character- 
ized by a strong clump effect with the increase of the average square distance that 
is slower than the Richardson law. These trajectories form structures, which are 
similar with fluid vortices and represent eddying regions. The size and the built- 
up time of the structures depend on the value of the initial potential. Trajectory 
structures appear with all sizes, but their characteristic formation time increases 
with the size. These structures or eddying regions are permanent in static stochas- 
tic potentials. The saturation time Tg represents the average time necessary for 
the formation of the structure. In time dependent potentials the structures with 
Ts > Tc are destroyed and the corresponding trajectories contribute to the diffusion 
process. These free trajectories have a continuously growing average displacement 
and dispersion. They have incoherent behavior and the clump effect is absent. The 
probability distribution functions for both types of trajectories are non-Gaussian. 

The average size of the structures S{K) in a time dependent potential is plotted 
in Figure 1. One can see that for i^T < 1 the structures are absent {S = 0) and that 
they appear for X > 1 and continuously grow as K increases. The dependence on 
X is a power low with the exponent dependent on the EC of the potential. The 
exponent is 0.19 for the Gaussian EG and 0.35 for a large EG that decays as 
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Figure 1: The average size of the trajectory structures for Gaussian EC (dashed hne) 
and for an EC that decays as (continuous hne). 

4.2 Anomalous diffusion regimes 

Test particle studies connected with experimental measurements of the statistical 
properties of the turbulence provide the transport coefficients with the condition 
that there is space-time scale separation between the fluctuations and the average 
quantities. Particle density advected by the stochastic E x B drift in turbulent 
plasmas leads in these conditions to a diffusion equation for the average density with 
the diffusion coefficient given by the asymptotic value of Eq. (4) . Recent numerical 
simulations [11] confirm a close agreement between the diffusion coefficient obtained 
from the density flux and the test particle diffusion coefficient. Experiment based 
studies of test particle transport permit to strongly simplify the complicated self- 
consistent problem of turbulence and to model the transport coefficients by means of 
test particle stochastic advection. The running diffusion coefficient D{t) is defined 
as the time derivative of the mean square displacement of test particles and is 
determined according to Eq. (4) as the time integral of the Lagrangian velocity 
correlation (LVC). Thus, test particle approach is based on the evaluation of the 
LVC for given EC of the fluctuating potential. 

The turbulent transport in magnetized plasmas is a strongly nonlinear process. 
It is characterized by the trapping of the trajectories, which determines a strong 
influence on the transport coefficient and on the statistical characteristics of the 
trajectories. The transport induced by the E x B stochastic drift in electrostatic 
turbulence [12] (including effects of collisions [13], average flows [14], motion along 
magnetic fleld [15], effect of magnetic shear [16]) and the transport in magnetic tur- 
bulence [17], [18] were studied in a series of papers using the decorrelation trajectory 
method. It was also shown that a direct transport (an average velocity) appears in 
turbulent magnetized plasmas due to the inhomogeneity of the magnetic fleld [19]- 
[21]. This statistical method was developed for the study of complex processes as 
the zonal flow generation [22], [23]. 

The results of all these studies are rather unexpected when the nonlinear effects 



are strong. The diffusion coefficients are completely different of those obtained in 
quasilinear conditions. A rich class of anomalous diffusion regimes is obtained for 
which the dependence on the parameters is completely different compared to the 
scaling obtained in quasilinear turbulence. All the components of particle motion 
(parallel motion, collisions, average flows, etc.) have strong influence on the diffusion 
coefficients in the non-linear regimes characterized by the presence of trajectory 
trapping. 

The reason for these anomalous transport regimes can be understood by analyz- 
ing the shape of the correlation of the Lagrangian velocity for particles moving by 
the E X B drift in a static potential [24] . In the absence of trapping, the typical LVC 
for a static field is a function that decay to zero in a time of the order Tfi = \c/V- 
This leads to Bohm type asymptotic diffusion coefficients Db = V'^^fi = VXc- Only 
a constant c is influenced by the EC of the stochastic field and the diffusion coeffi- 
cient is D — cDb for all EC's. In the case of the E x B drift, a completely different 
shape of the LVC is obtained for static potentials due to trajectory trapping. A 
typical example of the LVC is presented in Figure 2. This function decays to zero 
in a time of the order tji but at later times it becomes negative, it reaches a min- 
imum and then it decays to zero having a long, negative tail. The tail has power 
law decay with an exponent that depends on the EC of the potential [12]. The 
positive and negative parts compensate such that the integral of L{t), the running 
diffusion coefficient D{t), decays to zero. The transport in static potential is thus 
subdiffusive. The long time tail of the LVC shows that the stochastic trajectories in 
static potential have a long time memory. 

This stochastic process is unstable in the sense that any weak perturbation pro- 
duces a strong influence on the transport. A perturbation represents a decorrelation 
mechanism and its strength is characterized by a decorrelation time r^. The weak 
perturbations correspond to long decorrelation times, > Tfi. In the absence of 
trapping, such a weak perturbation does not produce a modification of the diffusion 
coefficient because the LVC is zero a,t t > Tfi. In the presence of trapping, which 
is characterized by long time LVC as in Figure 2, such perturbation influences the 
tail of the LVC and destroys the equilibrium between the positive and the negative 
parts. Consequently, the diffusion coefficient is a decreasing function of t^. It means 
that when the decorrelation mechanism becomes stronger (r^ decreases) the trans- 
port increases. This is a consequence of the fact that the long time LVC is negative. 
This behavior is completely different of that obtained in stochastic fields that do not 
produce trapping. In this case, the transport is stable to the weak perturbations. 
An influence of the decorrelation can appear only when the later is strong such that 
Td < Tfi and it determines the increase of the diffusion coefficient with the increase 
of Td- This inverse behavior appearing in the presence of trapping is determined by 
the fact that a stronger perturbation (with smaller Td) liberates a larger number of 
trajectories, which contribute to the diffusion. 
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Figure 2: Typical Lagrangian velocity correlation in static potential. 

The decorrelation can be produced for instance by the time variation of the 
stochastic potential, which produces the decay of both Eulerian and Lagrangian 
correlations after the correlation time Tc- The decorrelation time in this case is 
Tc and it is usually represented by a dimensionless parameter, the Kubo number 
defined by Eq. (3). The transport becomes diffusive with an asymptotic diffusion 
coefficient that scales as Dtr = cVXcK'^, with 7 in the interval [-1, 0] (trapping 
scaling [12]). The diffusion coefficient is a decreasing function of Tc in the nonlinear 
regime K > 1. 

For other types of perturbations, their interaction with the trapping process pro- 
duces more complicated nonlinear effects. For instance, particle collisions lead to 
the generation of a positive bump on the tail of the LVC [13] due to the property of 
the 2-dimensional Brownian motion of returning in the already visited places. Other 
decorrelation mechanisms appearing in plasmas are average component of the ve- 
locity like poloidal rotation [14] or the parallel motion that determines decorrelation 
when the potential has a finite correlation length along the confining magnetic field. 
The effects of an average component of the velocity are discussed in Section 5.1. in 
connection with drift turbulence. 

5 Test modes on drift turbulence 

Test particle trajectories are strongly related to plasma turbulence. The dynamics 
of the plasma basically results from the Vlasov-Maxwell system of equations, which 
represents the conservation laws for the distribution functions along particle trajec- 
tories. Studies of plasma turbulence based on trajectories were initiated by Dupree 
[8], [9] and developed especially in the years seventies (see the review paper [7] and 
references there in). These methods do not account for trajectory trapping and thus 
they apply to the quasilinear regime or to unmagnetized plasmas. A very important 
problem that has to be understood is the effect of the non-standard statistical char- 
acteristics of the test particle trajectories on the evolution of the instabilities and of 
turbulence in magnetized plasmas. 
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We extend the Lagrangian methods of the type of [9], [25], [26] to the nonhnear 
regime characterized by trapping. We study hnear modes on turbulent plasma with 
the statistical characteristics of the turbulence considered known. The dispersion 
relation for such test modes is determined as function of the characteristics of the 
turbulence. We consider the drift instability in slab geometry with constant mag- 
netic field. The combined effect of the parallel motion of electrons (non-adiabatic 
response) and finite Larmor radius of the ions destabilizes the drift waves. 

The gyrokinetic equations are not linearized around the unperturbed state as in 
the linear theory but around a turbulent state with known spectrum. The perturba- 
tions of the electron and ion distribution functions are obtained from the gyrokinetic 
equation using the characteristics method as integrals along test particle trajectories 
of the source terms determined by the average density gradient. 

The background turbulence produces two modifications of the equation for the 
linear modes. One consists in the stochastic E x B drift that appears in the trajec- 
tories and the other is the fluctuation of the diamagnetic velocity. Both effects are 
important for ions while the response of the electrons is approximately the same as 
in quiescent plasma. They depend on the parameters of the turbulence. 

5.1 The statistics of the characteristics 

The solution for the potential in the zero Larmor radius limit is 

0(x,z,i) = 0o(x- V*i,z), (9) 

where 0o is the initial condition and V^, is the diamagnetic velocity. This shows that 
the potential is not changed but displaced with the diamagnetic velocity. The finite 
Larmor radius effects consist in the modification of the amplitude and of the shape 
of the potential, but this appears on a much slower time scale. 

The ordering of the characteristic times for the drift turbulence is 

t\^t^^Tc^t\, (10) 

where r|, Ty are the parallel decorrelation times for electrons and ions (Ty'* = Ay/v^^* 

with Aj| the parallel correlation length and v^l^ the thermal velocity), r* = Ac/K is 
the characteristic time for the potential drift and Tc is the correlation time of the 
potential. The linear and nonlinear regimes are determined by the position of the 
time of flight in this ordering. The latter is much smaller than Tc and much larger 
than rjj. The statistical characteristics of the trajectories essentially depend on the 
ratio T^/Tfi. 

The quasilinear case corresponds to t^/t^i <^ 1 {V/V^ <^ 1), which means tur- 
bulence with the amplitude of the E x B drift smaller than the diamagnetic velocity. 
The motion of the potential produces in this case a fast decorrelation and trapping 
does not appear. The probability of displacements is Gaussian and the diffusion 
coefficient is Dqi = V'^t^. 

The nonlinear case corresponds to T^/rfi > 1 {V/V^, > 1). The motion of the 
potential is slow and trajectory structures produced by trapping exist in this case. 

The test particle motion in a drifting potential is obtained by a Galilean trans- 
formation from the motion produced by a stochastic E x B drift and an average 
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velocity Vd- This process was studied in [17]. It was shown that strips of opened 
contour hues of the effective potential + xVd appear due to an average velocity Vd 
and that the width of these strips increases with Vd until they completely eliminate 
the closed contour lines (for Vd> V). The Lagrangian correlation of the velocity in 
the presence of an average velocity Vd < V does not decay to zero as in Figure 2, 
but it has a positive asymptotic values at t — > oo. Consequently the transport along 
the average velocity is super diffusive in the static potential and diffusive with large 
diffusion coefficient (proportional with the average velocity) in the time dependent 
case. A part of the particles are trapped and the other move on the strips of opened 
contour hues of the effective potential. The invariance of the distribution of the 
Lagrangian velocity shows that the average velocity of the free particles Vj^ fulfils 
the condition 

nfrV}, = Vd, (11) 

and thus it is larger than the average velocity {Vj^ > Vd). ritr is the fraction of 
trapped trajectories and nfr is the fraction of free trajectories at a moment {ritr + 
Ufr = 1). 

This physical image leads, by changing the reference frame, to the following 
paradigm of the statistics of trajectories produced by the E x B drift in a moving 
potential. The trapped particles (structures) are advected by the moving potential 
while the other particles have an average motion in the opposite direction with a 
velocity Vfr such that 

rifrVfr + ritrV* = 0, (12) 

which is the equivalent of Eq. (11). This shows that there are particle flows in op- 
posite directions induced by the drifting potential if the amplitude of the stochastic 
E X B velocity is larger than the velocity of the potential. This determines a spliting 
of the probability of displacements in two parts: the probability of trapped and the 
probability of free particles. The first is a picked function that has constant width 
and moves with the velocity K- The second, is a Gaussian like function with an av- 
erage displacement {X2) — Vfrt — —V^t ritr/rifr- The probability of displacements 
at t < Tc is modeled by 

P(x, y, t) = ntrG(x, y - VJ; S^, Sy) + nfrG(x, y - Vfrt; + 2D J, Sy + 2Dyt) (13) 

where G{x,y; S^, Sy) is the 2-dimensional Gaussian distribution with dispersion 
Sx,Sy. We have considered for simphcity the distribution of trapped particles as 
a Gaussian function but with small (fixed) dispersion. The shape of this function 
does not change much these estimations. The free trajectories have dispersion that 
grows linearly in time. 

5.2 The growth rate of drift modes in turbulent plasma 

The average propagator of for a mode with frequency uj and wave number k = 
{ki,k2) is evaluated using the above results on trajectory statistics. It depends on 
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the size S{K) of the structures and on the fractions of trapped and free particles: 



/dr (exp (— ik • x(r))) exp (iuj(t — r)) 
-oo 



= iexp{-k'fS^) 



ntr _^ rifr 



u + kyV* u + kyVfr + ikfDi 



(14) 



where x(r) is the trajectory in the moving potential integrated backward in time 
with the condition x at time t. 

The solution of the dispersion relation is obtained as 



oo = kyV:" (15) 

(16) 



T.e// ^ V ro.^(n/r -ntr)+ Ifltr 
2 - To^ 



T^^-^[-\klS^^j (17) 

where Fq = exp{—b)Io{b), b = k\p\/2 and is the ion Larmor radius. The tensor 
Rij has the dimension of a length and is defined by 

Rji{T,t) = de' deM^i{\e\) (i9) 

J T J — OO 

where Mij is the Lagrangian correlation 

Mj,{\e' - e\) = {vj (x^(^'), ^, 0') d^vi (x^(^), z,e)), (20) 

and Vj is the E x B drift velocity. 

Several effects appear in the test modes characteristics due to the background 
turbulence. The spreading of ion trajectories produces the diffusion Di that influ- 
ences the growth rate (18) both in linear and nonlinear conditions. This term is 
similar to the result of Dupree, but the values of is influenced by trapping. 
Beside this, there are several influences that appear only in the nonlinear regime. 
The first is the factor given by Eq. (17), which is produced by the trajectory 
structures. It determines essentially the modification of the mode frequency. The 
flows of the ions induced by the drifting potential are represented by the fractions 
ritr and n/^-The tensor Rij is determined by the fluctuations of the diamagnetic 
velocity due to the background turbulence. We will analyze each of these processes 
in the next section. 



6 Trapping effects on the test modes 

The trajectory trapping has a complex influence on the mode. This can be under- 
stood by considering the evolution of the drift turbulence starting from a stochastic 
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potential with very small amplitude as it can be deduced from the growth rates of 
the test modes. 

The trajectories are Gaussian, there is no trapping in such potential and the 
only effect of the background turbulence is the diffusion of ion trajectories that 
produce resonance broadening. The well known results of drift modes in quasilinear 
turbulence are obtained 

. = *,v.^. , = ^(M±z^^iipW-*ji),. (21) 

Z — io \f^z\^Te z — io 

where = Dy = Dgi = V'^Xc/V^. This shows that the modes with large k are 
damped due to ion trajectory diffusion as the amplitude potential increases. The 
maximum of the spectrum is for u — kyVt,/2 and corresponds to k^pi^ ~ 1. 

When the nonlinear stage is attaint for V >V^^ the first effect is produced when 
the fraction of trapped trajectories is still small by the quasi-coherent component of 
ion motion. The structures of ion trajectories determine the T factor (17), which 
modifies the effective diamagnetic frequency (16) and the frequency uj. At this stage 
the flows can be neglected {ntr = 0, Ufr = 1) and Rji = in Eqs. (15)-(18), so that 
only the factor JF is important. It is interesting to note that this factors appears 
in Eqs. (15)-(18) as a multiple of Fq, although they come from different sources (JF 
from the propagator and Fq from gyro-average of the mode potential). This shows 
that the trapping or eddying motion has the same attenuation effect as the gyro- 
average. The maximum of the spectrum that appears for uj — k,yV^:/2 is obtained 
for smaller fcj^, of the order of the size of the trajectory structures k±S ~ 1. This 
means that the unstable range of the wave numbers is displaced toward small values. 
The maximum growth rate is not changed but displaced at values of the order 1/ S. 
Consequently, in this stage both the amplitude of the turbulence and its correlation 
length increase. 

At larger amplitude of the background potential, when the fraction or trapped 
ions becomes comparable with the fraction of free ions, the ion flows induced by 
the moving potential become important. These flows determine the increase of 
the effective diamagnetic velocity (16) toward the diamagnetic velocity and the 
modification of the growth rate of the drift modes. The latter decreases and for ritr — 
Ufr it is negative. The evolution of the amplitude becomes slower and eventually 
the growth rates vanishes and changes the sign. Thus, the flows of the ions induced 
by the moving potential produce the damping of the drift modes. 

The fluctuations of the diamagnetic velocity due to background turbulence de- 
termine a direct contribution to the growth rate (the tensor Rij). This term is zero 
for homogeneous and izotropic turbulence and strongly depends on the parameters 
of the anizotropy. The i = j = 1 component corresponds to zonal flows. Preliminary 
results show that it appears for trapped particles due to the anizotropy induced by 
ion flows with the moving potential. 



7 Summary and conclusions 

We have discussed the problem of stochastic advection of test particles by the E x B 
drift in turbulent plasmas. We have shown that trajectory trapping or eddying have 
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complex nonlinear effects on the statistical characteristics of the trajectories and on 
the transport. The nonlinear effects are very strong in the case of static potentials. 
The trajectories are non-Gaussian, there is statistical memory, coherence and they 
form structures. These properties persist if the system is weakly perturbed by time 
variation of the potential or by other components of the motion (collisions, poloidal 
rotation, parallel motion). The memory effect (long tail of the LVC) determines 
anomalous diffusion regimes. 

The process of trajectory trapping also influences the evolution of the turbulence. 
Recent results on test modes on turbulent plasmas are presented. They are based 
on a Lagrangian method that takes into account the trapping or eddying of the 
ions. The growth rate and the frequency of the drift modes on turbulent plasmas 
are estimated as function of the characteristics of the turbulence. The effects of 
the background turbulence appear in particle trajectories (characteristics of Vlasov 
equations) and in the fluctuations of the diamagnetic velocity produced by the den- 
sity fluctuations. We show that the nonlinear process of trapping, which determines 
non-standard statistical properties of trajectories, has a very strong and complex 
influence on the evolution of the turbulence. It appears when the amplitude of the 
E X B drift becomes larger than the diamagnetic velocity. 

A different physical perspective on the nonlinear evolution of drift waves is ob- 
tained. The main role is played by the trapping of the ions in the stochastic potential 
that moves with the diamagnetic velocity. We show that the moving potential deter- 
mines flows of the ions when the amplitude of the E x B velocity is larger than the 
diamagnetic velocity. A part of the ions are trapped and move with the potential 
while the other ions drift in the opposite direction. These opposite (zonal) flows 
compensate such that the average velocity is zero. The evolution of the turbulence 
toward large wave lengths (the inverse cascade) is determined by ion trapping, which 
averages the potential and determines a smaller effective diamagnetic velocity. The 
ion flows produced by the moving potential determine the decay of the growth rate 
and eventually the damping of the drift modes and generate zonal flows due to their 
nonlinear interaction with the fluctuations of the diamagnetic velocity. 
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